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The influence of repeated projective measurements on the dynamics of the state of a quantum 
system is studied as a function of the time lag r between successive measurements. In the limit 
of infinitely many measurements of the occupancy of a single state the total system approaches a 
uniform state. The asymptotic approach to this state is exponential in the case of finite Hilbert space 
dimension. The rate characterizing this approach undergoes a sharp transition from a monotonically 
increasing to an erratically varying function of the time between subsequent measurements. 

PACS numbers: 03.65.Ta, 03.65.Xp, 03.65.Aa, 05.30.-d 



I. INTRODUCTION 

Even though the theory of measurements provides an 
integral part of Quantum Theory and even though to this 
day it has passed any close scrutiny the measurement pro- 
cess remains a topic of controversial discussions. Accord- 
ing to von Neumann [T] the measurement of an observable 
yields one of the eigenvalues of the hermetian operator 
representing the measured observable and projects the 
state of the system onto the eigenspace of the resulting 
eigenvalue. There exist several attempts to explain this 
state reduction as resulting from the interaction of the 
system with a measurement apparatus whereby the sys- 
tem combined with the environment undergoes a unitary 
time evolution. Depending on the special structure of 
the measurement apparatus and its coupling to the sys- 
tem the reduced system state undergoes decoherence [21- 
[S] with respect to the eigen-basis of the operator to be 
measured. At the same time the state of the measure- 
ment apparatus adopts a particular "pointer state" . The 
latter process, also known as "einselection" , makes it pos- 
sible to read out the result of the measurement [frSj . 

As for other physical processes, any measurement re- 
sulting from an interaction of the system with the mea- 
surement apparatus will take a finite amount of time. 
This aspect is not taken into account in von Neumann's 
idealized picture where a projective measurement is con- 
sidered as an instantaneous process. It is therefore con- 
ceivable to consider arbitrarily frequent measurements 
which will lead to a total freezing of the dynamics, a phe- 
nomenon that is known as the Zeno effect [HHHj. When 
a finite, but still sufhciently short time r elapses between 
two measurements of the same observable the probability 
to consecutively find the system in the same eigenstate of 
the observable exponentially decreases with a rate that 
is proportional to the time r separating two measure- 
ments. Experimental verifications of the Zeno effect are 
reported in Refs. [151 HH]- Most recently, direct conse- 
quences of projective measurements have been discussed 
in different contexts such as thermodynamic control |17j , 
quantum tomography fl8l , and fluctuation theorems |19j . 



In this work we investigate the combined action of un- 
observed quantum motion and repeated projective mea- 
surements. In contrast to the standard treatment of the 
Zeno effect, according to which the system is followed 
only until it has left the initial state for the first time 
we keep track of all possible outcomes up to a total ob- 
servation time. In this way we take into account also 
realizations with intermediate deviations from and pos- 
sibly multiple revivals of the initial state. This approach 
differs from the standard treatment of the Zeno effect 
imposing selective measurements U, which disregards 
the future fate of the system after it has left the ini- 
tial state for the first time. We find for systems with 
finite-dimensional Hilbert space that asymptotically, in 
the limit of infinitely many measurements, all states of 
the system acquire equal probability. Hence, the proba- 
bility to retrieve the initial state is given by the inverse 
of the Hilbert space dimension. This result is universal 
and in particular does not depend on the time r between 
two consecutive measurements as long as r is different 
from zero. Exceptions from this universal behavior are 
found only for those singular values of the time between 
measurements and of the system parameters for which an 
invariant state exists of that part of the system that is 
not directly effected by the measurement. The approach 
to the stationary state eventually becomes exponential. 
The corresponding decay rate is a monotonically increas- 
ing function of r up to the time scale r* corresponding to 
the inverse of the largest energy difference of the system. 
For times larger than t* the decay rate shows a highly 
irregular behavior. 

This paper is organized as follows: Sec. II is devoted to 
the introduction of the system and the measurement of 
interest. In Sec. Ill, we investigate the time evolution of 
the density matrix under the repeated measurements and 
point out the universality of the stationary state of equal 
probability, followed by a more detailed discussion and 
proof in Sec. IV. We then focus on the survival prob- 
ability in Sec. V, and its decay behavior is detailed in 
Sec. VI. The paper ends with a summary and conclu- 
sions in Sec. VII. 



mwvwvv© 

BWVWVVO 

© 




FIG. 1: (Color online) Illustration of the system and the mea- 
surement protocol: The counter displays '1' when the particle 
occupies the left-most colored site, referred to as dot. When 
the particle resides at any chain site represented by one of 
the wells, the counter reports the event as '0'. Each time the 
one-handed clock has completed a period t a projective mea- 
surement is performed and its result displayed by the counter 
which remains idle in the meanwhile. 



II. THE SYSTEM 

We consider the dynamics of a single particle on a lat- 
tice consisting of a single site, referred to as dot and 
a one-dimensional chain with N sites acting as a reser- 
voir. A schematic view of this system at different stages 
during a measurement cycle is presented in Fig. 1. Any 
state of the particle can be represented in terms of lin- 
ear combinations of the discrete position states \x) where 
X = indicated by the coloring in Fig. [l] denotes the dot 
position and x = i with £ — 1 . . . N the chain positions, 
respectively. The total system is described by the Hamil- 
tonian 



N 



N-1 



H = Y^ e,\e}{£\ - 7c E (1^) (^ + 1| + 1^ + l)(^l) 



e=i 



(1) 



■7(|0)(1| + |1>(0|), 



where ei denotes the on-site energy of the site i relative 
to the dot-site energy eo = 0. Transitions between neigh- 
boring sites of the chain are governed by the parameter 
7c which will be used as energy unit, hence 7c = 1. Here 
we assume that only one end {£ = 1) of the chain is 
coupled to the dot. This model Hamiltonian has various 
applications. On one hand, it can be interpreted in terms 
of the Wigner-Weisskopf model POJ for the decay of the 
excited state of a two-level atom coupled to a finite num- 
ber of electromagnetic modes of an empty cavity. On 
the other hand, it describes a quantum dot coupled to a 
one-dimensional wire and, moreover, the dynamics of a 
cold atom in a one-dimensional optical lattice. Recently 
the decay of an unstable quantum state was discussed 
in terms of the Hamiltonian M with vanishing on-site 
energies e^ = 0, identical hopping rates 7c = 7 and in- 
finite N |5T]. In contrast to our present study effects of 



repeated measurements were not taken into account in 
this study. 

We consider the measurement of the occupancy of the 
dot-state that is represented by the projection operators 
onto the dot (d) and the chain (c), respectively. 



N 



H,-|0)(0|, Hc = 5]|f)(€| 



(2) 



e=i 



Note the completeness, n^ + Hc = 1, and the orthogonal- 
ity, nj = Hj, i = c,d, H^Hc = of the projection oper- 
ators. For chains with more than a single state, TV > 1, 
the observation of the dot-state occupancy corresponds 
to a selective measurement in the sense that the state 
reduction only removes dot-chain correlations but leaves 
the internal chain part of the state unchanged. 



III. DENSITY MATRIX 

We denote the density matrix of the total system im- 
mediately after the Af-th measurement by pj^j. The 
discrete time dynamics transforming the density matrix 
from the moment immediately after the M-th measure- 
ment to the moment immediately after the M -I- 1-st mea- 
surement, 



Pm+1 = MpM , 
is governed by the linear map Ai given by 

Mp= J2 nc«C/(T)p[/^(T)H, 



(3) 



(4) 



a^c,d 



It is composed of the unitary time evolution generated by 
U(t) = exp {—iHr/h) and a subsequent measurement- 
induced state reduction onto the subspaces of dot and 
chain. This reduction cancels any quantum correlation 
between dot and chain at each measurement by setting 
the nondiagonal matrix elements (0|p|£) with ^ = 1 . . . A^ 
to zero. Hence, all solutions of (Is]) are of the form 



Pai = PA/Hd + Xa/Hc 



(5) 



where pm is a real number, while xm is a reduced density 
matrix of the chain. Equation ([3| then decomposes into 
the following recursion relations for p and x: 

PM+i - ml(o|t/(T)|o)p + (o|c/(T)xAfC/t(T)|o) (6) 

XM+i = PAfHcC/(T)HrfC/t(T)Hc-hHcC/(r)xAfC/^(T)Hc, 

which can be iteratively solved. As the initial condition 
we choose Po = 1 and the null matrix Xo = 0- In passing 
we list the following properties of the map ^4 without 
proof: A^ is a trace-preserving, completely positive and 
contractive linear map |22j . These properties imply that 
the spectrum of Ai is contained in the unit disk of the 
complex plane. 
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FIG. 2: (Color online) (a) The tomography of the density 
matrix for a homogeneous chain (7 = 1, e^ = 0) with A'' = 9 
sites, displays the absolute values of the matrix elements at 
two times T. The unitary time evolution of the system is 
interrupted by measurements at equal intervals of duration 
T = 1. In the left panel at the comparatively short time 
T = 4 the diagonal elements are unevenly occupied and most 
of the non-diagonal elements have clearly non-zero values. 
The right panel illustrates the state reached at T = 50. Then 
the density matrix has become almost diagonal with equal 
weights, 0.1, and almost vanishingly small off-diagonal ele- 
ments in good agreement with the fixed point solution ([7|. 
(b) Density matrix evolution for a dot coupled to a chain 
with A'^ = 4 sites, 7 = 1, and random on-site energies taken 
from a uniform distribution within [—0.5,0.5]. The diagonal 
components Pdiag converge to 0.2, and the average magnitude 
of the off-diagonal components, poff defined in the text, decays 
exponentially in time (see the inset). 



We present tomographic views of the density matrix 
Pm in Fig. [21 a) for different times T = Mr with the time 
between measurements r = 1, as result of the M-fold ap- 
plication of the map M . Here and in the sequel we mea- 
sure time in units of l/(/i7c). After a few measurements, 
T — A, almost all elements of the density matrix are dif- 
ferent from zero, including most of the non-diagonal ones. 
The diagonal elements assume different values indicating 
a non-uniform population of the dot and chain sites. At 
T = 50, a stationary state with almost vanishing non- 
diagonal and identical diagonal elements has been ap- 
proached. These results represent the typical behavior as 
also exemplified by the panel (b) of Fig. [2] visualizing the 
dynamics of a system with on-site energies randomly cho- 
sen from a uniform distribution. As a measure of the off- 
diagonal elements, pij — {i\p\j), we determined the aver- 
age of their absolute values Pos = ^i^j \Pi,j\/Wi^ + ^)]: 
see the inset. Also in this case, we find the density matrix 




FIG. 3: (Color online) The diagonal elements of the density 
matrix in the position representation are encoded according 
to the scale given on the right as functions of position £ on 
the vertical axis and time T on the horizontal axis. Here, a 
dot is coupled to a chain with 50 sites, ee = and 7=1. For 
frequent measurements (r = 0.2 in the upper panel), the Zeno 
effect is visible by the persistent population in the dot (see 
also the enlargement), whereas in the lower panel for t — 5 
the dot population decays fast along the front ray that moves 
with the maximal group velocity 2 of the chain. 



approaching the uniform stationary state, as in Fig. K[a) 
for identical on-site energies. 

A more detailed account of the time-evolution of the 
diagonal elements is presented in Fig. [3] displaying the 
propagation of the population from the dot through the 
chain. In the upper panel with frequent measurements 
(t = 0.2), the transfer of the population from the dot to 
the chain is hampered as a result of the Zeno effect, see 
also the enlarged detail on the top of the figure display- 
ing the dot state and the first few chain sites. Behind 
a clearly visible front determined by the fastest signal 
propagation in the chain, a gradually decreasing excita- 
tion profile emerges. In contrast, for t — 5 the popula- 
tion behind a pronounced front ray decays fast on average 
with a superimposed periodic structure, which, at later 
times after the first reflections of the leading front, de- 
velops speckles. These speckles disperse only at much 
later times. It is also notable that independent of the 
time T the population dominantly propagates through 
the chain with the velocity 2. That can be understood 
from the fact that, apart from local effects at the border 
to the dot, the dynamics is dominated by the unitary 
time evolution in the chain. For a chain decoupled from 
the dot, the corresponding Hamiltonian M with ei — 
and 7 = 0, can be diagonalized to yield the dispersion re- 
lation, £{k) = — 2cosfc, where k — mr/N is the discrete 
momentum with n = 1,2,--- ,N. The group velocity 
is given by Vg = d£{k)/dk = 2sinfc, and its maximum 
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FIG. 4: (Color online) The decay of the dot-state sur- 
vival probability P = Pt/t in presence of a homogeneous 
chain with A^ = 9 sites is displayed as a function of the 
time T = Mr after M measurements for different values of 
r; the other parameters are the same as in Fig. [2] Panel 
(a) presents frequent measurements with short intervals, r = 
0.01, 0.02, 0.05, 0.1, 0.2, 0.5 from the upper to the lower curve. 
For these small r values the decay becomes faster with in- 
creasing r. Panel (b) presents the difference of P from the 
asymptotic state P* =0.1 for larger values of r. In this 
regime the decay rate does no longer monotonically increase 
with T. In all cases at large times T an average exponential 
decay sets in, which in some cases, is superimposed by an 
oscillatory behavior. 
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Poo = 



N +1 



(7) 



presents an invariant solution, where 1 is the unit matrix. 
In the Appendix |A] we prove that other solutions may 
only exist if there is an invariant chain state, i.e. a state 
with Mp = p and HdP = pTld = 0. 

This situation is exceptional; typically any chain state 
will evolve a dot component upon its unitary evolution. 
We note that the proof of the uniqueness of (IT]) holds 
true for general systems with finite-dimensional Hilbert 
space provided the measured observable possesses a non- 
degenerate eigenstate. The approach to the universal 
state (It]) therefore is not restricted to the special model 
(fTl). Most notably, the uniform state ^ corresponds to 
the thermodynamic equilibrium at infinite temperature. 
Apparently the repeated process of measurements "heats 
up" the system to an extreme degree. 



V. DOT POPULATION 

The relevant aspects of the dynamics in approaching 
this extreme state are displayed by the time evolution 
of the survival probability of the dot population which 
is experimentally easier to access than the whole density 
matrix. After M measurements, this survival probability 
is given by the dot matrix element pM = (0|/9m|0). Ac- 
cording to (Is]) and (Jl]), pj\/ does not only include the un- 
interrupted sojourn on the dot, but also allows for those 
cases in which the particle leaves the dot and reenters it 
at a later time. 

For a reservoir with a single level, that is, N = 1, 
analytic results are available. Since then lie = |l)(l|j the 
survival probability is simply given by the products of 
the 2x2 transition matrix T as 



PM = (T*Oii 



(8) 



with 



value is then Vg = 2. Hence it is the maximum group 
velocity that determines the front rays. 



Too Toi 
Tio Til 



(9) 



IV. EQUAL STATIONARY PROBABILITY 



The emergence of stationary states induced by re- 
peated measurements was also reported for special sys- 
tems [231 m] ) but up to now, neither conditions for their 
existence have been known, nor their universality has 
been noticed. In the following we address these ques- 
tions. The stationary solutions of the master equation 
([3]) are invariant under the action of the map Ai and 
therefore are solutions of Aipoo — Poo- One shows by 



where Tij = \{i\U{T)\j)\'^. The unitary time evolution 
operator U{t) of this two-level system can be easily 
solved to yield Too = 1 — (7/fi.ri)^ sin^(rir), where the 
oscillation frequency is given by fi = hr^ ^^"^ + 6^/4. 
Using Tqi = Tio resulting from the unitarity of U{t), 
and the probability conservation of the unitary evolu- 
tion, Too +2^10 = Til +T01 = 1, upon diagonalizing T we 
obtain 



PA/ = ^[i + (2roo-i) 



Ml 



(10) 



For the measurement with a very short time interval, 
Ht ^ 1, one finds Too ~ 1 + ©(r^), and accordingly. 
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FIG. 5: (Color online) The rate F = -In(!Ai|)/7^ deter- 
mined by the eigenvalue Ai that in absolute value is closest to 
1 is displayed as a function of the scaled time f = t/t* for dif- 
ferent values of the coupling constant 7 = 1, 0.75, 0.5, 0.25. 
Here r* denotes the shortest time scale of the unitary mo- 
tion. All curves start out from zero displaying a universal 
f^ growth behavior up to f ~ 0.4. Around f = 0.75 the F 
values decrease with decreasing 7. Above f = 1 the rate F 
presents an erratic behavior with large variations upon small 
parameter changes. 



Pm ~ 1 as long as M is not too large. This is a manifesta- 
tion of the quantum Zeno effect caused by frequent mea- 
surements that obstruct the time evolution of a quantum 
state. For larger values of r or after a large number of 
measurements, Mt'^ is no longer negligible; consequently 
Pm then deviates from unity and eventually converges to 
1/2, unless the measurement time interval accidentally 
exactly matches with the period, 27r/il, for which Tqq is 
one, see the expression for Too below Eq. Q. 

For reservoirs with more than one state numerical 
results are presented in Fig. p[a) where the survival 
probability P = pt/t is depicted as a function of the 
elapsed time T = Mr for relatively frequent measure- 
ments with accordingly small r. The smaller the time 
interval between the measurements is, the slower the de- 
cay of the survival probability becomes, in accordance 
with the Zeno effect. For instance, for the uppermost 
curve (r = 0.01), the survival probability stays close to 
unity even after 10^ measurements. On the other hand, 
for r = 0.5, the decay becomes more rapid and P quite 
soon converges towards its asymptotic value 1/10 which 
is in agreement with the survival probability resulting 
from for iV = 9. 

The panel (b) of Fig. [4] displays the deviation of the 
survival probability from its asymptotic value P* =0.1 
in the limit of large times. The decay rate sensitively 
depends on the time interval t although for larger values 
of T it no longer grows monotonically with increasing r 
values. In general, the approach to the asymptotic value 
crucially depends on the system details such as the size 
N and the coupling strength 7. 



VI. DECAY SPECTRUM 

The survival probability pM can be expressed in terms 
of the eigenvalues A„ and the eigenf unctions $„ and 
ifn of the linear map A4 and the corresponding dual 
map Ai^ , respectively, solving the eigenvalue problems 
M^n — A„$„ and M'^fn — ^nfn ■ The dual map is 
defined such that TiipM^ — Tr^ M^ip for ah matrices <I> 
and If. Using the bi-orthogonality of the eigenfunction, 
Tr(/3„$„, = dn,m and the initial condition pq = |0)(0| we 
obtain for the survival probability [25j 



PM ^Y.^n{0\^n\0){0\^n\0) 

n 

= l/(iV + l) + 5]Ar(0|^„|0)(0|<i>„|0), 



(11) 



where, in the second line, we split off the stationary part 
determined by Aq = 1, (^0 = 1 and $0 = l/(^ + 1)- 
The asymptotic decay towards the uniform stationary 
state is exponential with a decay rate given by the 
eigenvalue Ai whose absolute value is closest to unity. 
In practice though, the number of measurements to 
reach this asymptotic exponential decay regime may 
be extremely large due to the smallness of the coeffi- 
cient (0|(/Ji|0)(0|$i|0) and the presence of other eigen- 
values close to unity in absolute value. This may lead 
to extremely long transients with exponential decay at 
larger rates than those following from Ai, or even non- 
exponential decay laws. 

The eigenvalue Ai closest to unity depends in a reg- 
ular way on the length of the interval r as long as r is 
shorter than the fastest time scale t* « 27t/AE where 
AE is the width of the spectrum of the total Hamilto- 
nian. The rate defined as the negative logarithm of |Ai|, 
r = — ln(|Ai|)/7^, vanishes at r = and starts to grow 
proportionally to r^, in accordance with the Zeno effect. 
As can be seen from Fig. [5] the scaling of the negative 
logarithm of |Ai| by 7^ leads to a universal behavior of F 
at small r values. For values of r larger than r* the de- 
pendence of the rate F becomes very irregular and occa- 
sionally changes drastically from extremely small values 
corresponding to almost no decay of the survival proba- 
bility to rather large values giving rise to rapid decay. 

Finally we note without proof that a much simpler 
behavior emerges for an infinite chain with continuous 
energy spectrum and a spectral width which is much 
larger than any other energy scale. Under these condi- 
tions the dot population decays exponentially even in the 
absence of any measurement, i.e. |(0|C/(t)|0)| = e^*/"^" 
with To > 0. Then and only then a history in which 
the dot state is left at some time will never show a recur- 
rence of the dot population, or mathematically expressed 
(0|C/(T)ncC/(T)|0) = 0. Only in this case the traditional 
treatment of the Zeno effect yields the actual survival 
probability at a time t taking into account the full histo- 
ries up to this time. 



VII. SUMMARY AND REMARK 

We demonstrated that many measurements that are 
separated by finite periods of unitary time evolution lead 
to the universal infinite temperature state of equal prob- 
ability. We expect this behavior also to apply for more 
general systems with finite dimensional Hilbert space if 
they are subject to a large number of measurements of the 
occupancy of part of the system. This poses the question 
whether this result is a mere artifact of von Neumann's 
instantaneous reduction postulate, on which our analy- 
sis was based, or whether the repeated interaction with 
a measuring device drives the system into this extreme 
state. 

This work was supported by the National Research 
Grant funded by the Korean Government (NRF-2010- 
013-C00015). 



Appendix A: Proof of the uniqueness of the infinite 
temperature state 

Any normalized solution of the stationary master equa- 
tion 

p = Mp (Al) 

with the master operator M given by Eq. Q can be 
expressed as 



p^t/{N + l) + x, (A2) 

because the first term of the right hand 



where Tree 

side of Eq. (|A2|) already takes into account the full nor- 



malization of the density matrix p. We demonstrate that 
a nontrivial solution with x ^ Q and hence, other than 
infinite temperature solutions, may only exist if the uni- 
tary time evolution U{t) has invariant chain-states i.e. 
if 



y 



IiJJ{T)yU\T)Ii, 



(A3) 



has non-trivial solutions. Under typical conditions this 
will not be the case. The formal proof of Eq. (A3 1 as 



a necessary and sufficient condition is next presented. 
Since t/{N + 1) solves Eq. \A1\ one obtains 



: Mx . 



(A4) 



Hence, in analogy to Eq. (rol) x can be decomposed into 
a dot and a chain part reading 



qlid + z 



(A5) 



where g is a real number and z an operator on the chain 
sub-space, hence , sat isfying IVdZ = zlVd = 0, with Trz = 
—q. Using Eq. (A4) one obtains coupled equations for 
the unknowns d and z 

z = n,C/(r)2[/t(T)n, + qUMiT)UdU\T)Uc (A6) 
qlld = IldUiT)zU\T)Ud + qndU{T)UdUHT)Ud . (A7) 



We separately consider the two possible cases with q — 
and q ^ 0. 



(a) For q = Eqs. ( |A6| ) and ( [A7| ) simplify to read 

z = n^U{T)zUHT)U^ (A8) 

= UdUiT)zU\T)Ud. (A9) 



The first condition (A8| has the form of Eq. iA3l and 



hence has solutions only under exceptional conditions on 
the unitary chain dynamics and the time r. We note 
that Eq. (A9) is a consequence of Eq. (A8|. This follows 



by taking the trace over both sides of Eq. ( AS I and by 
replacing z by U{t)zW{t) under the trace on the left 
hand side, yielding 

Tr [U{T)zU'f{T)] ^ Tr [U^U (t) zU'' {t)Uc] (AlO) 



and consequently 



Tr[ndC/(r)zf/t(r)] -0 



(All) 



which is equivalent to Eq. ( |A9[ ) because 11^ is the projec- 
tion onto a one-dimensional subspace. 

(b) If g ^ one can divide both Eqs. ( |A6|A7[ ) by q 
and obtains for z ~ z/q two inhomogeneous equations 
reading 

z = n,[/(T)zf7t(r)n, + n,u{T)ndU^iT)ii, (A12) 
Ud = UdU{T)zuHT)nd + ndUiT)ndU\T)nd . (A13) 

By similar arguments as above one finds that the sec- 
ond equation follows from the first one. It is therefore 
sufficient to only consider the Eq. (A12). One finds by 
inspection as solution of Eq. (A12) z — He- This however 



does not present an admissible solution because it does 
not satisfy the normalization condition Trz = —q imply- 
ing Trz = -1 in contrast to Trz = Trllc = N. Whether 
there exist further solutions of the inhomogeneous equa- 
tion (A12| can be decided by means of the Fredholm 
alternative. It says that more than one solution of an 
inhomogeneous equation may only exist if the homoge- 
neous part of the equation possesses nontrivial solutions 
and if the inhomogeneity is orthogonal to all solutions 
of the dual homogeneous problem, see Sec. |VI[ One can 
show that in the present case the condition for the inho- 
mogeneity is automatically satisfied. The homogeneous 
part of Eq. (A12| again leads to a condition of the form 
of Eq. (A3). Assuming the existence of a nontrivial nor- 
malized solution u (Tru = 1) of ((A3| we obtain for the 



general stationary solution of Eq. ( Al I the expression 



P 



1 



l + N 



+ q] l-q{N + l)u, 



(A14) 



where u is a chain-state operator (JldU = ulld ~ 0) and 
q can assume arbitrary values in {—1/{N + 1), N). 

Hence in both cases (a) and (b) the uniqueness of the 
solution of the stationary master equation is guaranteed 
if and only if Eq. (A3 ) only has the trivial solution y — 0. 
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